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Abstract 

It is shown that a sufficiently strong external electric field causes a decrease in the transition 
temperature of ferroelectric, antiferroelectric, and metal-insulator transitions. The temperature 
dependence of the critical electric field suppressing a low-temperature phase is obtained. The 
relation between the low-temperature charge gap (the band gap) in the insulating phase and the 
metal-insulator transition temperature is also inferred. 
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In insulators, due to the presence of localized ionic charges, a strong external electric field 
produces mechanical stresses. It was shown recently [1] that a sufficiently high mechanical 
stress or tension causes a decrease in the melting temperature of a crystalline solid. Here 
we show that mechanical stresses produced by a sufficiently strong electric field in a low- 
temperature insulating phase lead to a decrease in the transition temperature of ferroelectric, 
antiferroelectric, and metal-insulator transitions. Such a decrease has been observed recently 
for the Verwey transition in magnetite (Fe304) from a high-temperature "bad metal" con- 
ducting phase to a low-temperature insulating phase [2]. A linear decrease in the Curie 
temperature caused by a strong external electric field is known for the antiferroelectric tran- 
sition in lead zirconatc (PbZrO^) [3]. In the case of metal-insulator transitions, the ratio 
of the low-temperature charge gap (the band gap) in the insulating phase to the transition 
temperature can be inferred from a general criterion of phase transitions in crystalline solids 
based on a concept of the critical number density of elementary excitations [4] . 

It was shown recently [1] that in a crystalline solid under a sufficiently strong mechanical 
stress a — F/S, applied to a sohd, where F is the applied force and S is the cross-section 
area of the solid in the plane perpendicular to the direction of the applied force, the energy 
Ey of the vacancy formation decreases in accordance with the formula 

^ E^Q - aa/no, (1) 

where Ey^ is the energy of the vacancy formation in absence of mechanical stresses, tiq ^ 1.1 x 
lO^^cm"^ and a ~ 18 are quantum constants, and a is the absolute value of the mechanical 
stress. A thermodynamic consideration based on the Clausius- Clapeyron equation gives 
the number density n„ of vacancies in a solid in the form [1] 

n„ = (Po/T) exp (-EjT) = (rioTo/T) exp {-EjT) , (2) 

where E^ is the energy of the vacancy formation, Pq — noTo is a constant, Tq can be put 
equal to the melting temperature of the solid at ambient pressure (in absence of mechanical 
stresses), and the Boltzmann constant kB is included in the definition of the temperature T. 

The melting of a crystalline solid occurs when the critical number density ric of vacancies 
is achieved, the critical number density of vacancies being [4] 
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Tie = riQexp (—a) . (3) 

The equations ([2]) and ([3]) give the relation between the energy of the vacancy forma- 
tion and the melting temperature Tm of a crystalline solid in the form 

= aT^. (4) 

A decrease in the energy of the vacancy formation and, whence, in the melting temper- 
ature is produced by the strong atomic relaxation in a crystalline solid under a sufficiently 
strong mechanical stress. 

This concept can be extended on other phase transitions, such as ferroelectric, ferromag- 
netic, and superconducting, if we consider the corresponding phase transition as the melting 
of some lattice [4]. In the case of the ferroelectric transition, it is the lattice of electrical 
dipoles; in the case of the ferromagnetic transition, it is the lattice of magnetic dipoles. 
Consider now the number density n of elementary excitations in a corresponding lattice, 
instead of the number density of vacancies in an ordinary lattice. An elementary excitation 
invokes the change of the sign of the electrical or magnetic dipole moment of the elementary 
lattice cell for ferroelectrics and ferromagnets, respectively. 

Similarly to the energy of the vacancy formation ([1]), the energy E of an elementary 
excitation depends on the absolute value aof the mechanical stress (in the region of strong 
mechanical stresses) as follows 

E = Eq - apcr/riQ, (5) 

where Eq is the energy of an elementary excitation at ambient pressure (in absence of 
mechanical stresses), and ap is a dimensionless constant dependent on the type of the phase 
transition [4]. The atomic relaxation constant ap is related to a dimensionless pressure 
derivative of the transition temperature Tc (in the region of high pressures) as follows 

k = riodTc/dP = —ap/a. (6) 

Here P is the hydrostatic pressure and the values of k and ap for ferroelectrics are 
k ^ —1/9 and ap ^ —ka ^ 2[4]. 
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It can be inferred from the experimentally measured pressure derivative of the metal- 
insulator transition temperature in BaVS^ [5] that k ^ —1/9 and ap ~ —ka ^ 2 for 
metal-insulator transitions too. Note that both ferroelectric and metal-insulator transitions 
are associated with a lattice distortion of a ferroelastic type (see, e.g., [6]). 

The number density n of elementary excitations is given by the formula analogous to the 
equation ([2]), 

n=inoTjT)expi-E/T), (7) 

where denotes here the transition temperature at ambient pressure (in absence of me- 
chanical stresses and external fields). 

The corresponding phase transition occurs when the critical number density Uc of ele- 
mentary excitations (given by the equation ([3])) is achieved [4], so that the energy E of an 
elementary excitation is related to the transition temperature T^., similarly to the equation 

iD, 

E ^ aT,. (8) 

In the case of metal-insulator transitions, the energy E of an elementary excitation is 
equal to the charge gap A^h (the band gap Eg). Optical reflectivity measurements in BaVS^ 
[5] show that the ratio of the charge gap Ach = Eg = E to the metal-insulator transition 
temperature Tmi is constant along the phase equilibrium curve in the temperature-pressure 
plane, in agreement with the relation ([H]). Since the charge gap is temperature dependent, 
it is a low-temperature value Ach (0)of the charge gap which should be compared with the 
value of the charge gap given by the equation ([8]). In BaVS^ at ambient pressure, the ratio 
of the charge gap to the metal-insulator transition temperature is A^h (0) /Tmi — 17 at 
10 K and Ach/Tmi = 12 at 60 K, as it can be determined from the low-frequency optical 
conductivity data reported in Ref.5. (The transition temperature Tmi at ambient pressure 
is Tmi ~ 70 K). 

In germanium Ge, the band gap is equal Eg (0) = O.TMeV, and the relation gives the 
metal-insulator transition temperature Tmi — 500K. 

Multiferroic bismuth ferrite (BiEeO^) exhibits a metal-insulator transition from the 
metallic cubic 7 phase to the insulating ferroelastic orthorhombic (3 phase at Tj\// ^ 1200/^ ~ 
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O.lOSeV^ [7]. The relation ([8]) gives a band gap Eg = aTui — l.S5eV. Optical absorption 
measurements show that the band gap of bismuth ferrite decreases slowly and linearly with 
temperature in the a and (3 phases from 2.5eVto 1.5eV [7]. (The rhombohedral a phase 
which is stable below 1100 K is ferroelectric and undergoes an antiferromagnetic transition 
at Tn ^ 6A0K). 

In cadmium chalcogenides, the band gap Eg is proportional to the melting temperature 
Tm or, equivalently, to the energy E^ of the vacancy formation, Eg = (3/4) E^: 



Material 


T K 


E^, eV 


Eg, eV 


Eg/Ey 


Crystal 
structure 


CdS 


2020 


3.13 


2.42 


0.77 


hexagonal 


CdSe 


1541 


2.40 


1.74 


0.73 


hexagonal 


CdTe 


1365 


2.12 


1.56 


0.74 


cubic 



Similar is the ratio Eg/E^ in diamond (0.8) and a — AI2O3 (0.7). In silicon (Si) and 
germanium (Ge), the band gap is Eg = 0.4i?„. 

We assume that a sufficiently strong electric field E applied to an insulating material 
produces a mechanical stress a which depends linearly on the electric field E as given by the 
formula 



The coefficient P will be estimated below. Substituting expression ([9]) in the equation ([5]) 
for the energy E of an elementary excitation we obtain 



Now the equations ([7]) and ([3]) give the phase equilibrium curve in the temperature-electric 
field plane in the form 



a = f3E. 



(9) 



E = Eq — ap(3E/nQ. 



(10) 



{Eo - apPE/rio) /T = a ^ Eq/T,. 



(11) 
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The last equation gives the temperature dependence of the critical electric field Ec sup- 
pressing a low-temperature phase as follows 

E,-E,(0)(1-T/T,), (12) 
where the value of the critical electric field at zero temperature is 

E, (0) ^ (noEo) I («p/?) . (13) 

A mechanical stress a produced by a sufficiently strong electric field E in an insulating 
material with ionic charges ^^e in the unit cell (e is the charge of an electron and i is the 
number of an ion) can be estimated as follows 

a = ZeE/al, (14) 
where Z is the total sum of absolute values of ionic charges in the unit cell (in units of 

^ 1/3 

the electron charge e), Z = \zi\, and ao = Uq = 0A5nm has an order of the lattice 
parameter a. Here it is assumed that an external electric field acts only on a charged layer 
at the surface of an insulating material (in the bulk the field is much weaker). 

In lead zirconate {PbZrO^), a linear temperature dependence of the critical electric field 
for the antiferroelectric transition with the Curie temperature Tc ~ 500K is indeed observed 
[3], the extrapolated value of the critical electric field at zero temperature being Ec (0) = 
800kV/cm. In this case Z = 12per formula. 

The equation f|T^ gives an estimation of the coefficient (3 in the form 

p ^ Ze/al (15) 
Substituting this expression in the equation (IT^ . we obtain 

Ec (0) ^ (a/apZ) {Tjea^) . (16) 

In the case of lead zirconate {PhZrO'i)T(. ^ 500K,Z = 12, and the equation (1161) gives 
the value E^ (0) = 750^1^/ cm which is close to the extrapolated experimental value of the 
zero-temperature critical electric field indicated above. 

An estimate (IT^ is valid only for sufficiently strong electric fields. In the case of lead 
zirconate, E ^ 20kV/cm] for weak electric fields, E < 20kV/cm, the Curie temperature Tc is 
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independent of the electric field E. In general case, the mechanical response in weak electric 
fields depends on the crystal symmetry [3]. The relation f|T4|) cannot apply to the metallic 
conducting phase above the metal-insulator transition. 

Polarization in thin films of BiFeO^ grown on SrTiO^ substrates [7] has an order of 

^ (3/ (47r) = {Z/A-k) (e/a^) ^ e/al = ^^p^Cjcm^, (17) 

where Z = 12for BiFeO^. The critical electric field given by the equation (fT6|) is much 
weaker than the effective polarization field (displacement) AnPr. This rules out a direct 
effect of an external electric field on the energy of an elementary excitation. 

In PbZrO^, an external electric field E ^ 20kV/cm produces a transition from the 
paraelectric cubic phase to the ferroelectric rhombohedral phase above the antiferroelectric 
Curie temperature Tc ~ 500K and a transition from the antiferroelectric orthorhombic 
phase to the ferroelectric rhombohedral phase below the antiferroelectric Curie temperature. 
Sufficiently strong electric fields suppress the ferroelectric phase as it follows from the above 
consideration. This conclusion is supported by the phase diagram of Pb {Zr, Ti) O3 system 
[3]. The Ti doping in PhZrO^ acts as an external electric field stabilizing the rhombohedral 
ferroelectric phase. However, a sufficiently high Ti doping suppresses this rhombohedral 
phase. 

In BaTiOs, a weak external electric field (E < IkV/cm) produces an increase in the 
ferroelectric Curie temperature [3], but stronger electric fields suppress a ferroelectric phase, 
in accordance with the above consideration. 

The critical electric field E^ (0) = 800kV/cm in PbZrO^ has an order of the breakdown 
field in thin BiFeO^ films [8]. The breakdown field can be estimated as follows. According to 
the equations ([1]) and ([9]), a sufficiently strong external electric field E produces a decrease in 
the energy E^j of the vacancy formation and, whence, the decrease in the melting temperature 
Tjn of an insulating material. The melting curve of an insulator in strong electric fields is 
determined by the equation similar to the equation ( ITTl) . 

(^,0 - aPE/no) /T = a^ E,^/T^. (18) 

The last equation gives the temperature dependence of the melting field (or the decom- 
position field) Em in the form 
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^ E„ (0) (1 - T/T„) , (19) 
where the zero-temperature value of the melting field is 

E„ (0) ^ rioE^o/aP ^ n^Tj^. (20) 
Substituting an estimate (fT5|) for the coefficient fi in the equation (!20l) . we obtain 

E„ (0) = (1/Z) (T„/eao) . (21) 

In the case of bismuth ferrite BiFeO^, the decomposition temperature is Td ^ O.llel^ 
[7] and Z = 12 per formula, so that the equation ( 12T1) gives E^ (0) = 200A; V^/cm. The 
breakdown field in thin BiFeO^ films is about Ej, = 700kV/ cm [8] . The electrical breakdown 
in BiFeO^ thin films is associated with a phase separation of BiFeO^ into magnetite Fe^O^. 
The coercive field in BiFeO^ thin films is also very high, about three times larger than in 
other ferroelectric thin films [8]. The breakdown field for Sm-doped bismuth titanate thin 
films is about E;, = 500kV/cm [9]. 

To summerize, we obtained the dependence of the energy of an elementary excitation 
on the electric field for ferroelectrics, antiferroelectrics, and paraelectric insulators in the 
region of strong electric fields. We obtained the temperature dependence of the critical 
electric field suppressing a low-temperature phase for ferroelectric, antiferroelectric, and 
metal-insulator transitions. An estimate of the breakdown field in insulating materials is 
given. The relation between the charge gap (the band gap) in the insulating phase and the 
metal-insulator transition temperature is also obtained. 
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